Energy-scales convergence for optimal and robust quantum transport in photosynthetic complexes 
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Underlying physical principles for the high efficiency of excitation energy transfer in light-harvesting com- 
plexes are not fully understood. Notably, the degree of robustness of these systems for transporting energy is not 
known considering their realistic interactions with vibrational and radiative environments within the surround- 
ing solvent and scaffold proteins. In this work, we employ an efficient technique to estimate energy transfer 
efficiency of such complex excitonic systems. We observe that the dynamics of the Fenna-Matthews-Olson 
(FMO) complex leads to optimal and robust energy transport due to a convergence of energy scales among all 
important internal and external parameters. In particular, we show that the FMO energy transfer efficiency is 
optimum and stable with respect to the relevant parameters of environmental interactions and Frenkel-exciton 
Hamiltonian including reorganization energy A, bath frequency cutoff 7, temperature T, bath spatial correla- 
tions, initial excitations, dissipation rate, trapping rate, disorders, and dipole moments orientations. We identify 
the ratio of XT/'yg as a single key parameter governing quantum transport efficiency, where g is the average 
excitonic energy gap. 



PACS numbers: 

I. INTRODUCTION 

Life on the earth has been solar-powered via the mechanism 
of photosynthesis for four billions years [1]. Photosynthetic 
antenna complexes have evolved to harvest the sun's energy 
and efficiently transport it to reaction centers where it is stored 
as biochemical energy. The first few steps in photosynthe- 
sis represent highly sophisticated energy capture and transfer 
processes, as in an efficient solar cell. Specialized pigments 
in the antenna complexes absorb energy from sunlight, creat- 
ing electron-hole pairs known as excitons. The excitons then 
travel to reaction centers where their energy is converted and 
stored as chemical energy. The key feature for the success 
of photosynthesis is the high efficiency of exciton transport - 
as high as 99% in certain bacterial systems (SJ. Thus, photo- 
synthetic complexes provide an excellent model for designing 
efficient artificial excitonic devices. Unfoitunately, the fun- 
damental structural and dynamical processes contributing to 
such efficient migration of excitons are not fully understood. 

During the last two decades, many spectroscopic obser- 
vations for extended excitation states over multiple chro- 
mophores within light-harvesting complexes have been re- 
ported |3 , 4|. These results support the quantum mechanically 
delocalized nature of excitations in spatial coordinate. A va- 
riety of theoretical models have been proposed to capture the 
nature of delocalized quantum coherence effects in excitonic 
energy transfer ll5 HT6l . These models employ various gen- 
eralizations of Forster energy transport ifTTl [181 and Redfield 
theory |fT9ll in different perturbative limits. Specific effects of 
quantum coherence in multichromophoric donor to acceptor 
energy transfer rates were also studied in details [15J . Recent 
advances in multidimensional electronic spectroscopy provide 
evidence that long-lived quantum dynamical coherence can 
exist in exciton basis |20-28|. These experiments reported 
oscillatory coherences that were observed to last around 500 
/s, which is on the same order as the transport time-scale. 



The observations were made at cryogenic and physiological 
temperatures for several light-harvesting complexes including 
the FMO complex of green sulfur bacteria ||20l |26l , reaction 
center (RC) of purple bacteria |21|, light-harvesting complex 
II of higher plants |,22|, conjugated polymers ll24ll . and ma- 
rine algae ||251 . It was specifically suggested by Engel and 
Fleming that quantum dynamical interference effects might 
explain exceptionally high energy transfer efficiency of light- 
harvesting systems EOl . This conjecture has lead to vigorous 
theoretical efforts to try to understand the role of coherence in 
excitonic transport ||29H45 1. Despite considerable progress in 
experiment and theory, it is still unknown how quantum co- 
herence can be preserved in such unusually warm, complex 
and wet conditions. Moreover, it is not yet fully understood if 
quantum interference effects, in either spatial or energy coor- 
dinates, are contributing to efficiency of these systems. More 
importantly, it is not known how robust and optimal are theses 
complexes with respect to variations in the coherent system 
evolution and/or environmental parameters. 

Reference ||29l demonstrated that an effective collabora- 
tion between coherent quantum evolution and environmental 
fluctuations could enhance photosynthetic energy transfer ef- 
ficiency (ETE). This work was based on a quantum trajec- 
tory picture in site basis 1291 , within the Born-Markov and 
secular approximations which guarantee complete positivity 
of the excitonic dynamics B6]| . However, due to the per- 
turbative approximation used in this model, it was impos- 
sible to explore the optimality of energy transfer efficiency 
in the relevant regimes of intermediate system-bath coupling 
strength. Moreover, due to secular approximation, the coher- 
ence and population transfers are essentially decoupled. In or- 
der to avoid such problems subsequent studies relied on non- 
peiturbative Haken-Strobl method |6| to be able to explore 
a wider range of environmental interactions in the context 
of a (pure-dephasing) classical white-noise model at infinite 
temperature. These model studies illustrate optimal ranges of 
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dephasing-assisted excitation transfer ||3T1 -|331. In the optimal 
regime, an appropriate level of environmental fluctuations can 
wash out the quantum localization effects at the equilibrium 
state, when they are not too strong to lead to quantum Zeno 
effect |31|. These models, however, are by construction in- 
adequate to capture the role of quantum interplay of system 
evolution with non-equilibrium dynamics of bath within real- 
istic non-perturbative and non-Markovian regimes ||35|. 

Recently Ishizaki and Fleming developed a general ap- 
proach for studying excitation energy transfer in multichro- 
mophoric systems based on a hierarchy of coupled master 
equations |36|. Using this Hierarchy Equation of Motion 
(HEOM), first introduced by Kubo and Tanimura 1471, quan- 
tum coherent beating in the FMO protein at physiological 
temperature was investigated theoretically |37|. Although, 
HEOM can be in principle applied to all the different regimes 
of environmental interactions as a general benchmark, the re- 
quires computational resources significantly increase with the 
size of the system, the bath correlation time-scale, and within 
low temperature regimes. Thus, a number of alternative or 
complementary techniques for simulation of open quantum 
dynamics have been recently proposed |48-54|. We have re- 
cently shown that the second-order time-convolution (TC2) 
master equation can be used to efficiently estimate ETE in 
large complex excitonic systems interacting with bosonic en- 
vironments in the intermediate regimes |55|. This method is 
based on the Cao's earlier work on the generalized Bloch- 
Redfield master equations BSl . In Ref. Il55l we examined 
the TC2 master equation reliability beyond extreme Marko- 
vian and perturbative limits. 

In this paper, we observe that the FMO pigment-protein 
complex is optimal and robust with respect to all the estimated 
parameters of system and environmental interactions, within 
the TC2 master equation evolution for either of Lorentzian or 
Ohmic spectral densities of the environment. The paper con- 
tains the following main results: 

• We show that efficient and fault-tolerant energy transport 
occurs when the FMO internal parameters are in tune with en- 
vironmental parameters leading to a collaborative interplay of 
the coherent free Hamiltonian evolution and incoherent effects 
due to environment. 

• We comprehensively study the effects of reorganization en- 
ergy A, bath frequency cutoff 7, temperature T, positive and 
negative spatial correlations, loss, and trapping mechanisms 
on the energy transfer efficiency (ETE) landscapes. We show 
that ETE is optimal for the physiologically estimated value of 
these parameters, and is robust over a wide range of variations 
in these parameters. 

• We demonstrate that a convergence of time/energy scales for 
the relevant internal end environmental parameters of FMO 
complex facilitates efficient energy transport. For an exci- 
tonic system with average energy gap g, we observe that a 
single effective parameter A = \T/^g governs the behavior 
of quantum transport dynamics. Small values A ^ 1 give rise 
to weak localization and low efficiency. Intermediate values 



Awl correspond to the optimal energy transfer. Large values 
A 1 give rise to strong localization and low efficiency. 

• We observe that positive/negative spatial correlations es- 
sentially renormalize the reorganization energy to effective 
lower/higher values. The positive bath correlations can sig- 
nificantly enhance ETE at the regime of large A by inducing 
symmetries in the effective phonon-exciton Hamiltonian pro- 
tecting the transport against strong dynamical disorder In the 
intermediate values of system-bath coupling, the spatial cor- 
relations can enhance the robustness of transport efficiency. 
However, at the very small values of reorganization energy, 
they have an adversarial effect by diminishing useful but very 
weak bath fluctuations. 

• To explore the dynamical role of coherent system evolution, 
we examine whether or not the spatial structure of the FMO 
complex plays any important role in its performance. We find 
that the FMO compact structure is instrumental for support- 
ing high efficiency. Moreover, we directly demonstrate that 
the FMO performance is robust to variations in its geometri- 
cal parameters considering its rich internal parameter space. 
Specifically, the efficiency of FMO is robust to variations in 
the orientations of dipole moments and site energies due to its 
compact structure. 

• Finally, we discuss possibility of a so-called Goldilocks 
principle in the quantum regime that could explain the con- 
vergence of time-scales in FMO complex and other light- 
harvesting systems, with respect to internal parameters and 
environmental interactions. The Goldilocks principle for 
complex systems postulates a 'just right' level of complex- 
ity - too little complexity compromises function, while too 
much complexity compromises robustness. Convergence of 
timescales increases complexity by allowing different quan- 
tum processes to interact strongly with each other These in- 
teractions have the potential either to reduce or to increase 
efficiency and robustness. For the FMO complex and other 
simplified small-size light-harvesting systems, we show that 
the convergence of timescales effectively acts to increase both 
efficiency and robustness. Such a fundamental principle could 
also serve a key guiding rule for designing optimal materials 
to achieve exceptional quantum transport performance in re- 
alistic environments. 

The organization of the paper is as follows: section II re- 
views the TC2 master equation that we employ to efficiently 
simulate energy transport in multichromophoric systems in 
low excitation regimes [55]. In the section III, the efficiency 
and sensitivity of the FMO is studied with respect to the en- 
vironmental parameters including reorganization energy, bath 
cutoff frequency, temperature, bath spatial correlations, and 
trapping. In section IV, we consider the dependence of ETE 
on variations in temporal and geometrical factors of the trap- 
ping mechanism for FMO. The role of initial excitonic states 
is investigated in the section V. In the section VI, we study the 
robustness of ETE in presence of disorder in the FMO internal 
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FIG. 1: The disordered structure of the Fenna-Matthews-Olson 
(FMO) complex: It is a trimer consisting of tiiree identical monomers 
each formed from seven BChls embedded in a scaffold protein. The 
FMO complex acts as an energy transfer channel in green sulphur 
bacteria guiding excitons from the light-harvesting antenna complex, 
in the proximity of BChls 1 and 6, to the reaction center which is in 
the proximity of BChls 3 and 4. 

Structure parameterized by site energies, inter-chlorophyl dis- 
tances, and dipole moment orientations. Some complemen- 
tary materials for FMO complex are presented in the appen- 
dices on the FMO structural data and ETE in the presence of 
an Ohmic bath. 



II. THEORETICAL MODEL OF THE FMO COMPLEX 

To explore optimality and robustness of natural or artificial 
light-harvesting complexes we compute the energy transfer ef- 
ficiency landscape as a function of the independent degrees 
of freedom of the pigment-protein complex including thermal 
and vibrational environments. This imposes two fundamen- 
tal obstacles to developing analytical and computational ap- 
proaches that can be employed in this context. The first prob- 
lem is the fact that these systems are complex and contain a 
significant number of important physical variables. These in- 
dependent parameters and their dynamical interplay have to 
be fully explored as they could play a major role in the energy 
transfer dynamics within these complexes. The second prob- 
lem is that the common theoretical models, including Forster, 
Haken-Strobl, Redfield or Lindblad theory, are based on sim- 
plifying perturbative and/or Markovian assumptions that are 
inadequate by construction |35|. This second problem exists 
even if the estimated parameters of the system or system-bath 
couplings are rather weak and when bath has little or no mem- 
ory. The difficulty is that for studying the optimality or robust- 
ness of these complexes one has to fully explore both topology 
and geometry of energy transfer efficiency landscape outside 
of the known values of a given pigment-protein complex ||56J . 
The combination of these two problems indicates that such 
study will be practically intractable with current state of the 
art simulation techniques such as HEOM L37J . which, while 



accurate, are computationally highly expensive. Recently, we 
have derived the well-known time nonlocal master equation 
TC2 without making the usual weak system-bath coupling as- 
sumption. In this section, we summarize the main relevant 
definitions and equations of that approach. For more techni- 
cal details we refer the readers to Ref. 1551 . 

The dynamics of a multichromophoric system interacting 
with surrounding scaffold protein and solvent can be under- 
stood by starting from a general time evolution formulation 
of open quantum systems. The total system-bath Hamiltonian 
can be expressed as 

Htotal = Hs + Hph + Hs-ph (1) 

where 

Hs = Y.^j\]){j\ + J,k\j){k\, 
Hs-ph = ^SjBj. 



The phonon bath is modeled as a set of harmonic oscilla- 
tors [57'|. Here \j) denotes an excitation state in a chro- 
mophore spatially located at site j. The diagonal site energies 
are denoted by e^s that include reorganization energy shifts 
\j — huj^d'j due to interactions with a phonon bath; 
dj^^ is the dimensionless displacement of the (j, ^)th phonon 
mode from its equilibrium configuration. The strengths of 
dipole-dipole interactions between chromophores in different 
sites are represented by Jjk- The operators Si = and 
Bj = — f^^^dj^^Qj^c are system and bath operators. Here, 
we assume that each site is linearly interacting with a sepa- 
rate phonon bath. The overall dynamics of the system can be 
expressed by a time-nonlocal master equation, e.g., TC2, t55l 
as: 

^p(t) = Csp{t) + Ce-hp{t) (2) 
- Et^^' p ^^(* - t')e^''^'-''^S,pit')dt' - h.c] 

where Cj (t — ti) = {Bj {t)Bj (ti)) represent the bath correla- 
tion function, and the Liouvillian superoperators £5, Cph and 
^s-ph associated to Hs, Hph and Hs-ph respectively. The 

term Ce-h = - ^Ls{b')0'l' ■} ~ rtrap{\trap){trap\, .} 
captures two different competing electron-hole pair recombi- 
nation processes that determine the energy transfer efficiency 
of light harvesting complexes. The first process is the loss 
due to dissipation to the environment at each site that happens 
within the time-scale of 1 ns. This inevitable adverse envi- 
ronmental effect guarantees that the energy transfer efficiency 
has a value less than one. The second process is the desired 
recombination process due to successful trapping. 

The ETE is defined as accumulated probability of exciton 
being successfully trapped: 
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Hamiltonian parameters are given in appendix A, as functions 
of chromophoric distances, dipole moment angles and site en- 
ergies. The trapping rate is treated as a free parameter in sec- 
tion IV. For most of this manuscript we assume an excitation 
initially localized at site 1. Similar results are obtained for 
other initializations such as localized state at site 6. We study 
the overall, best, and worst case scenarios of dependency on 
the initial states in the section V. In the next section, we use 
Eq. (j3]l to demonstrate efficiency and robustness of the FMO 
complex with respect to variations in reorganization energy, 
temporal correlations, temperature and spatial correlations. 
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FIG. 2: The Energy transfer efficiency (ETE) of the Fenna- 
Matthews-Olson complex versus reorganization energy A (as a mea- 
sure of decoherence strength) and bath cutoff frequency 7 (as a mea- 
sure of non-Markovian character of the bath f65\ 1661 ). The ex- 
perimentally estimated values of T = 298 °K, A = 35 cm~^, 
"f — 50— 166 cm^^, r-(^,^p — 1 ps and r^J'^^ = Ins reside at an op- 
timal and robust neighborhood of ETE. The FMO complex would act 
sub-optimally in the regimes of large 7 and very small A, or large A 
and very small 7. In the intermediate values of A/7 the phenomenon 
of environment-assisted quantum transport takes place. A top view 
of this plot in Fig. [5] indicates A/7 is the parameter that governs 
efficiency at a fixed temperature. 
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r] = 2rtrap / {trap\p{t)\trap)dt (3) 
Jo 

The above performance function is biologically relevant 
and has been extensively used for a variety of light-harvesting 
complexes I9ll29ll55l. 

In this work, we mostly concentrate on the dynamics of the 
FMO pigment-protein complex of the bacterium Chlorobium 
tepidum [63 1 as a prototype for larger light-harvesting antenna 
complexes. Fig. [T] The structure of this pigment-protein com- 
plex was the first of light-harvesting complexes to be deter- 
mined by x-ray crystallography |62|. The FMO structure con- 
sists of a trimer, formed by three identical monomers, each in- 
cludes a closely packed assembly of seven BChl-a molecules. 
The FMO trimer guides excitation energy transfer between 
the chlorosome, the light-harvesting antennae of green sulfur 
bacteria, and the membrane-embedded type I reaction center 
(RC). In the following sections for our numerical simulations, 
unless specified otherwise, the environmental parameters for 
the FMO complex are chosen according to the experimentally 
or theoretically estimated values of reorganization energy 35 
CTO~^, temperature 298 °K, trapping rate of 1 ps, exciton life- 
time of 1 ns, and zero spatial correlations. In the present sim- 
ulations, we adopt the bath cutoff frequency of 50 cm~^ from 
the studies in Refs. |4, 72|. However we have repeated our 
simulations for another reported estimation of 7 = 166 cm^^, 
Ref |59|, and we observed no significant difference in the be- 
havior of ETE landscape. The diagonal and off-diagonal free 



III. OPTIMALITY AND ROBUSTNESS WITH RESPECT 
TO ENVIRONMENTAL PARAMETERS 

We employ the time-nonlocal master equation presented 
above to efficiently estimate the energy transfer efficiency 
landscape as a function of various independent system and 
environmental degrees of freedom over a wide range of val- 
ues. This efficient simulation allows us to examine compre- 
hensively all relevant regimes of the multiparameter space 
for finding possible high efficient and robust neighborhoods. 
Only after such exhaustive study can one quantify the per- 
formance of any particular natural photosynthetic complexes. 
Moreover, such studies shed light on the maximum capabili- 
ties that can be achieved for optimal material design to engi- 
neer and characterize fault-tolerant artificial light-harvesting 
systems ll60ll6Tl within a given rich system-environmental pa- 
rameter space. The quantum efficiency of photosynthetic en- 
ergy conversion can also be estimated by measuring the quan- 
tum requirements of ATP formation |2|. Thus, the overall ro- 
bustness and optimality of a pigment-protein complex can also 
be experimentally explored, verified, or calibrated by varying 
the tunable parameters in the laboratory. For example, this 
can be achieved by changing ambient temperature and using 
diverse solvents with different dielectric properties [61 1. 

In addition to our landscape study of ETE, we also ex- 
amine the degree of optimality and robustness of the energy 
transfer by employing the the Euclidean norm of the gradi- 
ent and Hessian matrix of the ETE function. The Euclidean 
norm of the ETE gradient at any parameters values pi and 
P2, ||V77(pi,p2)||2, quantifies the degree of optimality. The 
gradient measure reveals the degree of local optimality in a 
surface manifold. Careful inspections of the room tempera- 
ture plots for the ETE function versus various pairs of rel- 
evant parameters show a convex or concave manifold, thus 
gradient as a measure of local optimality suffices to measure 
global optimality. To examine the robustness, we compute the 
Hessian matrix norm I (77(^1, p2)) 1 12 = \\[d'^r]/dpidpj]\\2 
= 1,2) as the total measure of local curvature of the 
manifold. A smaller value of this norm corresponds to a flat- 
ter surface, thus a more robust process. We use a five-point 
stencil method to compute derivatives numerically. 
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FIG. 3: Top view of ETE landscape illustrated in Fig. [2] indicating 
that the ratio of the reorganization energy over bath frequency cutoff 
can be considered as the parameter that governs the energy transfer 
efficiency at a fixed temperature. As we inspect this plot in an angu- 
lar coordinate from the vertical axis, 7, toward the horizontal axis, 
A, we can distinguish different regions of the ETE landscape that are 
separated by straight lines A/7. At very small decoherence rate, the 
FMO complex experiences weak localization due to static disorder. 
As we increase this ratio, an optimal region of ETE emerges that is 
induced by an appropriate level of interplay between environmental 
fluctuations and coherent evolution. At higher levels of this param- 
eter, ETE drops significantly due to strong localization induced by 
dynamical disorder. 



A. Reorganization energy and non-Marliovianity 

We first explore the variation of the FMO energy transfer 
efficiency versus reorganization energy and bath cutoff fre- 
quency using a Lorentzian spectral density, see Fig. |2] The 
reorganization energy. A, is proportional to the squared value 
of the system-bath couplings, quantifying the decoherence 
strength. The bath cutoff frequency is the inverse of the bath 
coherence time-scale that captures the non-Markovian nature 
of the environment. That is, the non-Markovianity measure, 
defined as the information flow from the system to the phonon 
bath by Breuer et. al. Il65l . increases exponentially with de- 
creasing bath cutoff frequency |66 |. The ETE function has 
not been illustrated in Fig. [2] for bath cutoff frequency val- 
ues less than j — 5 cm^^ , since according to our analysis in 
Ref. 1 55 1 the simulation errors of TC2 master equation may 
become significant in such highly non-Markovian regimes. It 
should be noted that in the regime of 7 > /3~^ we apply low 
temperature corrections to the bath correlation function as ex- 
plained in the next section. 

The optimality and robustness of ETE for the FMO protein 
complex at the experimentally estimated values of A 35 
cm^^ and 7 = 50 — 166 cm^^ are evident in Fig. |2] An 
independent study on the optimality of ETE versus reorgani- 
zation energy has been recently reported in Ref. |67|. It can 
be observed in Fig. [2] that the non-Markovianity of the bath 
can slightly increase ETE in the regimes of weak system-bath 
coupling. However, such slow bath behavior can significantly 



decrease ETE when the system interacts strongly with it. The 
main question is how one can understand this phenomenon for 
all non-Markovian and Markovian regimes in the context of 
Environment-Assisted Quantum Transport (ENAQT), which 
was first investigated using simplified formalisms of the Lind- 
blad (weak coupling and Markovian assumptions) ll29ll and 
Haken-Strobl (infinite temperature, pure-dephasing assump- 
tion) ||30l. 

In this work, we report an underlying theory for 
environment-assisted quantum transport. The landscape in 
Fig. |2]shows a remarkable interplay of reorganization energies 
A and bath frequency cutoff 7. The ETE takes values below 
unity if the FMO operates at the limit of very small A and large 
7. On the other limit, the FMO efficiency drops significantly 
when operating at large A and very small 7. As we argue be- 
low these two regimes can be understood as manifestations of 
weak- and strong quantum localization, respectively. 

A note on terminology: we use the terms weak and strong 
localization in analogue to the corresponding effects in bulk 
solid-state systems |68|. However, in contrast to solid state 
systems, here we are dealing with finite-size systems of few 
chromophores. In addition, it is important to note that because 
of the long-range nature of the dipolar force, the bulk versions 
of these chromophoric systems would not exhibit exponen- 
tially localized states. What is really being investigated here, 
then is a kind of transient localization which would eventu- 
ally go away due to the long-range nature of the dipolar force 
and the small system size. If the destruction of transient local- 
ization takes longer than the exciton lifetime, however, then 
transient localization is just as adversarial to quantum trans- 
port efficiency as full localization. In this paper, when we re- 
fer to weak and strong localization, the reader should keep in 
mind that we are actually discussing weak and strong transient 
localization. 

In order to fully explore the regimes of weak and strong 
localization and the various intermediate transitional regions, 
we illustrate a top view of Fig. |2]in Fig. [3] This plot im- 
mediately reveals that the regions of distinct energy transfer 
efficiencies at room temperature are governed by a parameter 
of the form A/7. For small and for large values of this param- 
eter, the efficiency is low. The efficiency reaches its maximum 
for intermediate values of this parameter In both strong and 
weak quantum localization limits the excitation will be spa- 
tially trapped in the regions typically far from reaction center 
and eventually dissipating to bath due to adversarial electron- 
hole recombination processes which occur on three order of 
magnitude slower time-scale. In the intermediate regime, the 
right amount of the interplay of quantum coherence and envi- 
ronmental fluctuations can facilitate an optimal energy trans- 
port in a robust fashion by various physical mechanisms in- 
cluding minimizing site energy mismatches, washing out po- 
tential destructive quantum interference effects, and enhanc- 
ing the energy funneling by providing an appropriate vibra- 
tional energy sink. 

The dependence of efficiency on A/7 can be understood as 
follows. When the FMO complex interacts weakly with its 
environment the exciton migration is essentially dictated by 
the site energies and inter-pigment couplings. The spatial lo- 
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cations. Fig. [T] and dipole moments of BChls do not show 
any obvious regular pattern. In the absence of environmental 
interactions, destructive interference effects due to such ran- 
dom configurations can cause weak localization over a few 
chromophores away from the trapping site. This phenomenon 
prevents successful delivery of the initial exciton to the re- 
action center and therefore leads to below unity efficiency of 
about 83% at A = 0. Weak localization is amplified by low- 
ering the temperature, a behavior observed in the simulations 
presented in the next section. By increasing the system-bath 
coupling strength, the adversarial interference effects of exci- 
tonic pathways are diminished in a fashion similar to obser- 
vations reported in Refs. 1291 [30l [32l . However, here it can 
be observed in Fig. |3]that for larger bath cutoff frequencies, a 
larger reorganization energy is required for ENAQT to occur 

We can understand these observations in terms of ENAQT 
by noting that the effective decoherence rate is given by A/7 
in the perturbative limit at a fixed temperature ||46| . This ra- 
tio becomes smaller as we raise the 7 therefore a stronger 
coupling A is needed to guarantee the level of decoherence 
strength required for ENAQT. The overlap of a delocalized 
exciton wave function with the trap enables an almost com- 
plete, 98% quantum transport in the optimal range of A/7 at 
the ambient temperature. As we increase the reorganization 
energy and bath coherence time-scale, the ETE starts to drop. 
In this regime, excitonic wave function again experiences lo- 
calization as the environmental fluctuations change their role 
to play a strong adversarial effect on the quantum transport 
essentially as a source of strong dynamical disorder. The ETE 
landscape in Fig. [3] clearly has level sets that exhibit linear 
relationship with A and 7. Indeed the ratio A/7, known as the 
Kubo number, is the parameter that governs Anderson local- 
ization transition in stochastic classical modeling of environ- 
mental interactions. ||68147T1 . In the fixed high-temperature 
limit of ETE illustrated in Fig. |3] one can observe that A/7 is 
a determining parameter for transport efficiency in the regions 
beyond optimal ENAQT area. In the next section, we go be- 
yond the Kubo number, by directly investigating the tempera- 
ture dependent energy transfer dynamics, leading to a general 
governing parameter as Ar/7. 

A more quantitative study of the degree of optimality and 
robustness of the energy transfer as functions of system-bath 
coupling strength and bath memory is illustrated in Fig |4] 
The optimality is defined as Euclidean norm of the ETE func- 
tion gradient ||V7y(A,7)||2 to locate the local maxima in the 
ETE landscape in Fig. |4] The robustness is defined by 
||i/(7y(A,7))||2 to measure local curvature of the manifold. 
Note that the ETE gradient and Hessian matrix norms are in- 
dicated in a logarithmic scale, thus the global optimal point 
with zero derivative can not be explicitly highlighted in this 
representation. The experimentally estimated values for FMO 
are illustrated as black dots in each graph clearly located in an 
optimal and robust region. One remarkable feature is the fact 
that environmental parameters of FMO have almost the min- 
imal reorganization energy and bath cutoff frequency among 
all the regions with simultaneous optimal and robust proper- 
ties. This might be explained by the overall tendency in nature 
to minimize the amount of required work, that is, facilitating 



logdIVETEIlj) 




100 200 300 400 500 
Reorganization energy (cm 



100 200 300 400 500 
Reorganization energy (cm 



FIG. 4: This figure is a complement to Fig. |2] (Left) The degree 
of ETE optimality is quantified at different values of A and 7 by the 
gradient matrix norm of the ETE function. The dark blue points rep- 
resent near optimal values. (Right) The degree of ETE robustness is 
quantified by the Euclidian norm of Hessian of the ETE. The dark 
blue points represent near robust points. The estimated FMO en- 
vironmental values of A = 35cm^^ and 7 = 50 cm~^, marked 
by black dots, are located on the corner of both robust and opti- 
mal region. This implies the FMO environmental conditions provide 
the minimal work required for facilitating optimal and robust energy 
transport by maintaining a rather weakly-coupled small scaffold pro- 
tein. 



an optimal and robust environmental platform for the FMO 
energy transport by preserving a rather small size scaffold pro- 
tein that is weakly coupled. 

The theoretical and experimental modeling of the environ- 
ment surrounding the FMO complex suggests that the spectral 
density of the phonon bath modes can be expressed as a some 
of a few Lorentzian terms |57|. The spectral density function 
J{ijj) determines the time correlation functions as: 



{mm) = 



1 



dujj{uj)- ' (4) 

1 — exp(— pfiw) 



where a Lorentzian spectral function has the form J{lS) = 
2Aa;/ (w^ +7^)- For simplicity in this work we have consid- 
ered a single Lorentzian term with amplitude A = 35 cm~^ 
and cutoff frequency 7 = 50 cm~^^. However these val- 
ues were actually extracted by fitting the experimentally mea- 
sured absorption spectrum using an Ohmic spectral density 
J(a;) = A(w/7) cxp(— w/7) [72]. Thus, we also examine 
ETE versus variations of A and 7 in such a model depicted in 
Fig. 18 see appendix B. We note that for an Ohmic model a 



very similar behavior of the optimality and robustness of en- 
ergy transport can be observed as those in Fig. [2] however, 
the estimated ETE is lower than those obtained by exploiting 
a Lorentzian model. Specifically for FMO the ETE value are 
92.3% and 96.7% for Ohmic and Lorentzian spectral density 
respectively. This suggests that an Ohmic bath can not cap- 
ture the near unity efficiency of FMO complex, and favors the 
Lorentzian model as a more accurate description of the bath 
spectral density. 
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FIG. 5: Energy transfer efficiency manifold as a function of reorga- 
nization energy and temperature. The temperature range is from 35 
°K to 350° AT. At tfie regimes of large system-bath coupling strength 
the ETE drops faster by increasing the temperature due to strong lo- 
calization induced by dynamical disorder. There is a narrow region 
of reorganization energies between 10 and 50 cm^^ that ETE is rel- 
atively robust with respect to variations in temperature. The FMO 
value of A = 35 cm~^ is located in this region indicating a signifi- 
cant insensitivity to temperature variations. At very low reorganiza- 
tion energy and low temperatures weak localization can be observed 
due to static disorder. 



B. Reorganization energy and Temperature 

Here, we directly examine the temperature dependence of 
the effective parameter A/7 discussed in the pervious section. 
Generally light-harvesting complexes in bacterial, marine al- 
gae, and higher plants operate at various temperatures and re- 
organization energies. For example, the FMO complex of the 
green sulfur bacteria can operate at the bottom of the ocean 
at a depth of hundreds of meters, and also in hot springs with 
diverse temperature variations. Moreover, 2D electronic spec- 
troscopy of these systems have been performed at different 
cryogenic and room temperatures. Thus, we need to explore 
the efficiency and sensitivity of energy transport as a function 
of reorganization energy and temperature. 

We note that in the low temperature limit, (3 > 7^^, 
the Lorentzian bath correlation function {B{t)B{t'))ph is no 
longer a single term A(2//3 — 17)6^'*'*^*^* ^ and should be cor- 
rected by a sum of exponentially decaying terms 

where Vk — 2TTk/f3h are bosonic Matsubara frequencies. In 
practice a truncation of the above infinite series is needed for 
numerical simulations. The higher levels of truncation are dic- 
tated by lower temperature limits for the systems under study. 
To guarantee an accurate estimation of the bath correlation 




100 200 300 400 500 
Reargaization energy (cm 

FIG. 6: The top view of the ETE landscape presented in Fig. |5]as 
a function of inverse temperature and reorganization energy. In the 
high temperature regime (/3~^ > 7), the three regions of weak lo- 
calization, ENAQT and strong localization can be distinguished by 
the parameter AT given the free Hamiltonian of FMO and 7 = 50 
cm~^. These results, combined with Fig. [3] suggest that the pa- 
rameter Ar/7 governs the shape of the overall FMO energy transfer 
landscape. 



function for the calculation of energy transfer landscape, in 
the relevant low temperatures, we consider the first one hun- 
dred Matsubara frequencies in the above summation. It should 
be noted that the mentioned low temperature corrections have 
a minor numerical cost for calculating ETE by using TC2 
Q, since we just need to compute sum of the Laplace trans- 
forms of the above exponential terms |55|. This is another 
advantage of using the TC2 in contrast to the more accurate 
approaches such as HEOM for which every correction term 
iyj.e^"'''*"* ) necessitates considering N extra auxiliary vari- 
ables in the simulation, where N is the number of the sites, 
therefore significantly increasing the computational cost of 
the simulations. 

We show the behavior of ETE as a function of reorgani- 
zation energy A and temperature T in Fig. |5] It can be ob- 
served that at various possible conditions for the FMO pro- 
tein of T = 280 to 350 °X and A = 35 cm~\ the ETE has 
an optimal value and resides in a robust region of the energy 
transfer landscape. We note that for system-bath strength cor- 
responding to the FMO's environment, the energy transport is 
very robust to variations of temperature and just slightly de- 
creases in the low temperature limit. By increasing the reorga- 
nization energy, the temperature dependence of ETE becomes 
more pronounced. As we move toward the classical regime 
of the dynamics at high reorganization energy and high tem- 
perature ETE drops significantly. This implies the necessity 
of quantum effects for highly efficient and robust excitonic 
energy transfer It is remarkable that at the relevant physio- 
logical temperatures, the ETE is robust only within the range 
of A = 35 CTO^^, and it becomes unstable as we approach re- 
organization energy of over 100 cm~^. On the other hand, at 
the very low temperatures, enhancement of A would improve 
the ETE and brings it to a saturated high level. The gradient 



8 



and Hessian norms as functions of reorganization energy and 
temperature are illustrated in Fig. |7] At the relevant FMO 
operating temperatures, optimum and robust energy transport 
can occur simultaneously only within a small regions of A 
between 30 to 35 cm^^, which coincide with the estimated 
values of reorganization energy for the FMO. We note that 
there are certain regions of higher robustness at higher reor- 
ganization energy that are in principle available, but these re- 
gions imply a significantly lower operating temperature for 
the FMO operation and they have suboptimal ETE in compar- 
ison with actual FMO environmental parameters at the room 
temperature. 

As shown in the previous section, the parameter A/7 
governs the shape of the ETE landscape at a fixed high- 
temperature limit. However, in the perturbative limit the de- 
coherence rate can be expressed by XTf-f which captures the 
suboptimal ETE in the weak localization region. Now, we in- 
vestigate if the ETE behavior in all regimes can be globally 
captured by the parameter AT/7, for the given FMO Hamil- 
tonian. Specifically, we need to verify if one can predict the 
optimal noise-assisted transport region as well as ETE sup- 
pression levels at the strong localization regions by a sin- 
gle parameter XTf-f. To examine the validity of this the- 
ory, we study ETE as a function of the reorganization energy 
and the inverse temperature for the fixed 7 = 50 cm^^, see 
Fig. |6] Similar to the plot of efficiency as a function of A, 
7, Fig. [3] here also the efficiency landscape is divided by 
lines of approximately AT (with some deviations from lin- 
earity in the high A and low temperature regime). It can be 
observed from the Fig. |6]that for small AT weak localiza- 
tion is dominant. At the intermediate AT values environment- 
assisted energy transport occurs. As we move towards larger 
system-bath interactions and higher temperatures, strong dy- 
namical disorder diminishes the coherence and the exciton mi- 
gration can be fully described by an incoherent hopping pro- 
cess, since the wave function is essentially localized over spa- 
tial sites. At this regime, the effect of high temperature can be 
understood from the dynamics of BChls energy fluctuations 
which is described by the symmetrized correlation function 
S{t) = l{{B{t),B{0)})ph (Bj = B, for any BChl j). The 
function can be extracted experimentally by three-pulse pho- 
ton echo peak shift measurement. For a Lorentzian density, 
S{t) — 2hX/ f3e~''* , and fixed 7, the temperature and the re- 
organization energy determines the amplitude of the site en- 
ergy fluctuations. Our simulation in Fig. |5] demonstrates that 
a high efficient energy transfer can be achieved at a moder- 
ate site energy fluctuations away from both weak and strong 
locaUzation limits. 

Combining different regimes of three important environ- 
mental parameters. A, 7, and T, the effective decoherence 
strength AT/7 emerges as the parameter that governs the en- 
ergy transfer efficiency landscape. For the FMO free Hamil- 
tonian energy gaps, by increasing the single parameter AT/7 
from small to intermediate, and from intermediate to large val- 
ues, one can describe the transition from weak localization 
to ENAQT, and from ENAQT to strong localization. More 
generally, when the effective decoherence rate AT/7 is either 
much smaller or much larger than the typical energy splitting g 
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FIG. 7: This figure is a complement to Fig. |5] Left: the degree 
of ETE optimality is quantified at different values of A and T by 
the norm of the ETE gradient. The dark blue points represent near 
optimal values. Right: the degree of ETE robustness is quantified 
by the Euclidian norm of the ETE Hessian. The dark blue points 
represent near robust points. We note that within the range of pos- 
sible FMO operating temperatures (e.,g between T = 2SQ° K to 
T = 350° if) simultaneous optimal and robust energy transport can 
only be achieved for A values around 30 to 35 ann~^, that is equiv- 
alent to the estimated reorganization energy for the FMO complex 
that is marked by a black dot at room temperature T = 298° A". 



between delocalized energy eigenstates, then transport is sup- 
pressed. Thus, in order to predict the general patterns of quan- 
tum transport in generic light-harvesting systems, we should 
compare the relative strength of AT/7 (with dimension of en- 
ergy) to the average excitonic energy gap of the free Hamilto- 
nian that can be quantified as 5 = iv^ll^ ~ Tr{H)l/N\\^ 
for an N level system, where the nuclear norm is de- 

fined as the summation of X's singular values. This follows 
from the fact that the ETE of a system with a rescaled Hamil- 
tonian aH equals the ETE of a system with Hamiltonian H 
for which other environmental energy/time scales are renor- 
malized by a factor 1/a. Here we give a simple proof for this 
fact: 

The efficiency defined in Eq. Q can be simply calculated 
as 77 = 2rtrap{tfCLp\p{s = Q)\trap) where ^ denotes the 
Laplace transform. The operator p(s = 0) can be found by 
transforming the dynamical equation (|2]): 



p(0) - Csp[Q) + C,-hp{Q) 



(6) 



Y,[S,,^C{-Cs)Sfp{Q)~h.c.]. 



This yields the ETE 

?/ = -2rtrap{trap\}Cp{0)\trap) 

where 

1 p., 



(7) 



/C = (£5 + Ce-h - Y.^Sj, j^C{~Cs)Sj - h.c.])-^ (8) 
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FIG. 8: (a) ETE as a function of reorganization energy and positive bath spatial correlations that correspond to in-phase fluctuations of two 
spatially separate bath oscillators. At the large values of reorganization energy, the positive correlations induce certain symmetries in the 
effective phonon-exciton Hamiltonian significantly enhancing ETE, by protecting the system against strong dynamical disorder similar to 
decoherence free subspaces fl9]. In the intermediate values of A, the ENAQT become more robust due to spatial correlations. In large 
correlation length, the optimal ETE is essentially expanded over a much wider regime, one order of magnitude larger than the optimal ETE 
region in the absence of any spatial correlations. At smaller values of A positive spatial correlations actually hinder the exciton transport by 
decoupling useful but very weak fluctuations. A similar linear relationship is observed here between A and positive spatial correlations as those 
with 7 and T. (b) ETE as a function of reorganization energy and negative spatial correlations that correspond to out of phase fluctuations of 
two spatially separate bath oscillators. Negative correlations reduce ETE at all values of A. 



The Laplace transform of the correlation function (|5]l 
can be found as C{-Cs) = n(^^icot{^) - i) + 

Er=i Considering these explicit ex- 

pressions, it is easy to see that the system-bath parame- 
ters {ai/, A, 7, T, rtrap, r;oss} yields the same ETE (jvl) as 
the parameter set {H, A/a, 7/a, T/a, rtrap/a, rioss/ajTWe 
should mention that similar results can be obtained for the hi- 
erarchy equation of motion. 

Overall, we introduce the dimensionless parameter as A = 
AT/75, as the parameter that governs the energy transfer effi- 
ciency landscape. When we approach A « 1, the decoherence 
rate is tuned to give the maximum transport rate. For A ^ 1 
and A ^ 1 we move toward weak and strong-type of tran- 
sient localizations respectively. To examine the temperature 
and also spectral density independency of the observed pat- 
terns in Figs. |3]and[6j we present the landscape top view of 
ETE for a) Lorentzian density at 17° K and b) Ohmic density 
at 298° in appendix B. In both cases a similar pattern can 
be observed. 

In the subsequent sections, we fully explore the interplay of 
reorganization energy with a variety of Frenkel-exciton free 
Hamiltonians for systems consist of up to 20 chromophores 
with modest site energetic disorders. Our numerical studies 
suggest that in the optimal transport regime the relevant en- 
ergy scale of the free Hamiltonian, g, can be essentially ex- 
pressed as the average dipole-dipole interactions strength that 
is |/Lt|^ n, where is |^|^ is the average dipole-moment strength 
and n = is the density of N chromophores homogenously 
embedded in a system with size D. We will address the sig- 
nificance of the convergence of energy scales in the broader 
context of complex quantum systems in the discussion sec- 



tion. 



C. Reorganization energy and bath spatial correlations 

Recently there have been significant interest on the po- 
tential role of bath spatial correlations as the underlying 
physical principle leading to the experimental observations 
of long-lived quantum coherence in biomolecular systems 
II20I l2n l26l . These correlations have also been the center of 
attention and debate in the recent theoretical studies of quan- 
tum effects in photosynthetic complexes 11301 1741476 1 as a po- 
tential positive feature of environmental interactions. The ba- 
sic intuition is that if the environmental fluctuations are cor- 
related in space, e.g., in the site basis, then they will lead 
to fairly global quantum phase modulations in exciton basis, 
similar in spirit to decoherence-free subspaces that are stud- 
ied in details in quantum information science f79l. Indeed, 
as demonstrated earlier in Ref. |30|, using a Lindblad master 
equation, positive spatial coiTelation can enhance the under- 
lying contributions of quantum coherence to ETE according 
to two different measures based on Green's function methods. 
However, the overall ETE remains relatively unchanged for 
the FMO estimated values, as the contribution of relaxation to 
ETE drops with a similar rate Ref. Il30l . Here, we are inter- 
ested in the potential role of quantum correlations in optimal- 
ity and/or flexibility of energy transport rather than enhancing 
the time-scale of quantum coherence beating. 

The TC2 can be generalized to include the effect of both 
spatial and temporal correlations between protein environ- 
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FIG. 9: The degree of sensitivity of ETE for 10* uniformly sampled pure and mixed initial exciton density matrices for different values of 
reorganization energy: (a) The error bars indicate the ETE standard deviation. At the FMO value of A = 35cm~^ the standard deviation of 
ETE has a negligible value of about 0.1%. (b) A worst-case scenario of FMO energy transport sensitivity to initial exciton states. Here, the 
error bars indicate the maximum and minimum values of ETE achieved over the sample of 10* randomly chosen initial states. At the FMO 
value of A = 35cm~^ the variation in ETE is less than 1% 



ments: 

- Csp{t) + C,-hp{t) (9) 

ot 

- Et^'' ^ f " t'V'^'-''^s''p(t')dt' - h.c] 

j.k 

The function Cjj is just the autocorrelation function of the 
jth environment. The cross correlation between fluctua- 
tions of sites j and k is given by the function Cj^k{t ^ 
t') = {Bj{t)Bk{t')). In our study we assume Cj,k to be 
of Lorentzian Aj.fe7a;/(7^ + w^) form with the cut-off fre- 
quency 7, similar to the autocorrelation spectral density, and 
the strength Xj^k that decays exponentially with the distance 
between Bchls j and k, dj,k- A correlation length Rcor 
determines the strength of the spatial coiTelations, Xj ^ = 
X exp{—dj,k / Rcor) ■ We distinguish between the two differ- 
ent cases of positive and negative correlations, corresponding 
to A > and A < respectively. The sign of the correlations 
is imposed by the position and orientation of the Bchls with 
respect to their surrounding protein |77 1. 

The effects of positive and negative spatial correlations 
on the energy transfer efficiency for various degrees of the 
system-bath couplings are presented in Fig. [8] One main ob- 
servation is that both positive or negative spatial correlated en- 
vironmental fluctuations play an insignificant role in the over- 
all ETE only at the regions of reorganization energy very close 
to A = 35cm^^. However, as we slightly increase the system- 
bath coupling strength, we see two very distinct behaviors of 
the ETE as a function of long-range spatial characteristics of 
the bath fluctuations, depending on whether these variations 
being correlated or anti-coiTelated. For rather strong reorga- 
nization energy, positive spatial correlations can enhance the 
ETE by about 30%. By contrast, negative spatial correlation 
can reduce the ETE by the same amount for reorganization 
energy of above A = 400cm~^. 



A careful inspection of Fig. [8] (a) shows a linear relation- 
ship between reorganization energy and positive spatial corre- 
lations similar to those with 7 and T, as discussed in the pre- 
vious section. For very small reorganization energies, when 
quantum localization due to static disorders is present, the 
positive spatial correlations can slightly reduce ETE by de- 
coupling the fragile but positive role of environmental fluc- 
tuations. On the other extreme regime, positive spatial cor- 
relations can play a significant positive role by inducing an 
effective symmetry in the system-bath interactions, leading to 
partial immunity with respect to the strong dynamical disor- 
ders at large reorganization energies. Remarkably, in the in- 
termediate regime, the induced symmetries can substantially 
enhance the robustness of the environment-assisted transport. 
We can observe that the width of ENAQT region (in reorgani- 
zation energy axes) is enhanced from 20cm^^ for zero corre- 
lation length to 200cm^^ for 100 A correlation length. At this 
wide regimes, the system is protected with respect to the ad- 
versarial effects of dynamical disorder preventing the strong 
localization effects. Therefore, the spatial correlations should 
be also incorporated in the governing parameter AT/7. This 
can be achieved by renormalizing the reorganization energy 
A to an effective lower/higher values, when positive/negative 
correlations exist. 

In order to quantify the role of spatial coiTelations in en- 
hancing the quantum coherence time-scale of light-harvesting 
complexes, one has to partition the time-nonlocal master 
equation, Eq. |2] similar in spirit to those studies in Ref. 
Il30ll38l . In the other sections of this work we assume spatial 
correlations to be negligible for the FMO complex consistent 
with the recent simulations in Ref. |78|. Indeed, the atomistic 
simulation of the FMO pigment-protein-solvent dynamics at 
room temperature presented in Ref |78| reveals insignificant 
coiTelations in the site energy fluctuations, suggesting that the 
uncorrected bath approximations are reasonably valid. Next, 
we explore the role of initial conditions. 
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FIG. 10: The ETE manifold on the parameter space of reorganization 
energy and trapping time-scale. It can be seen that the FMO complex 
reaches its maximum functionality at trapping rates of about 0.5 — 5 
ps~^. The tunnel-shape ETE landscape can be understood by noting 
that at low trapping rates the transport efficiency is diminished by the 
recombination process. At high trapping rates the exciton transfer is 
suppressed via quantum Zeno effect as the trapping process corre- 
sponds to very strong and continuous measurement of the system. 



IV. ENERGY TRANSPORT SENSITIVITY ON THE 
INITIAL EXCITATIONS 

The exciton migration pathways and time-scales have been 
studied in detail for a variety of light-harvesting complexes 
using various perturbative techniques including Forster mod- 
els for studying LHI and II of purple bacteria |j9l and Lind- 
blad models for simulating the dynamics of the FMO protein 
of green sulphur bacteria [!29j. Nevertheless, the role of ini- 
tial conditions in the overall energy transfer efficiency of pho- 
tosynthetic complexes is to a large extent unknown. It was 
recently shown that the initial quantum coherence could influ- 
ence the energy transfer efficiency in LHI of purple bacteria 
assuming no interaction with the phonon bath |34|. The de- 
pendency on initial localized excitation at BChls 1 and 6 were 
also examined for the FMO complex using Lindblad, Haken- 
Strobl, and HEOM models |29, 31, 37 1. However, the sensi- 
tivity of ETE with respect to generic initial pure and mixed 
states taken from a large ensemble in the single-excitation 
manifold has not previously been explored. 

Here, we first examine the average sensitivity of ETE with 
respect to randomly chosen initial states for various reorgani- 
zation energies. To this end, for each value of reorganization 
energy, we sample over 10"' (pure or mixed) density matrices 
from a uniform distribution in the space of all 7 x 7 trace one 
positive matrices. In Fig. [9] (a) the average values of ETE is 
plotted with an error bar representing the variance of ETE in 
our random sampling. Note that at the optimal ETEs, corre- 
sponding to the value of reorganization energy of the FMO 
complex, the dependency of the variances on initial states is 
very small - less than 0.1%. However, the ETE fluctuations 
can grow by an order of magnitude for larger or smaller values 
of A. We also investigate the best and worst possible random 
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FIG. 11: This figure is a complement to Fig. [TO](Left). The degree 
of ETE optimality is quantified at different values of A and r^^\p 
by the norm of the ETE gradient. The dark blue points represent 
near optimal values. (Right) The degree of ETE robustness is quan- 
tified by the Euclidian norm of the ETE Hessian. The blue points 
represent near robust points. The FMO achieves its maximum effi- 
ciency at r^Jjp — 0.5ps. Note that for larger reorganization energies 
the trapping rate has be increased to achieve optimal ETE. However, 
this competition does not exist at small and intermediate system-bath 
coupling strength, where A is on the order of off-diagonal elements 
of the FMO free Hamiltonian leading to environment-assisted en- 
ergy transport. This suggests that a general convergence of time- 
scale might be required to obtain global efficient and robust transport 
(see discussion section). Although certain regions of highly robust 
ETE can be observed in figure (b), but those cases are suboptimal by 
considering their corresponding values in figure (a). 



initial single excitonic states in the Hilbert space of FMO. In 
Fig.|9](b), we illustrate this extreme possible deviations by er- 
ror bars on the average ETE at any given value of A. Note that 
ETE is very robust, varying about 1% with respect to different 
initial excitations at the optimal area of ETE landscape. How- 
ever, this robustness diminishes substantially at the regimes 
of large reorganization energy. Next, we study the ETE land- 
scape as a function of trapping and dissipation rates. 



V. TEMPORAL AND GEOMETRICAL EFFECTS OF THE 
TRAPPING MECHANISM 

Basic structural information on the FMO-RC complex has 
been obtained via linear dichroism spectra and electron mi- 
croscopy lISOl . These studies indicate that the symmetry axis 
of the trimer is normal to the membrane containing the re- 
action center. The electron microscopy resolution is generally 
not sufficient to distinguish the top and the bottom of the FMO 
trimer nor the distance between FMO-RC. Thus, in principle 
either pairs of BChls 1 and 6 or BChls 3 and 4 are the pig- 
ments that connect the FMO complex to RC. However, it is 
widely believed, due to efficient energy funneling toward RC, 
that BChl 1 and 6 are the linkers to antenna baseplate, and 3 
and 4 should serve as target regions within the neghibourhood 
of RC complex [|16J . This hypothesis has been recently veri- 
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FIG. 12: The ETE as a function of reorganization energy for the 
initially maximally mixed state. In each plot one of the 7 BChls is 
considered to be connected to the reaction center. BChls 3 and 4 
acting as the exciton transfer bridge yield the highest efficiency for 
almost all values of the reorganization energy. This confirms the 
experimental evidence that the FMO spatial orientation is such that 
the BChls 3 and 4 are located near the RC. 




FIG. 13: ETE versus dissipation (loss) time-scale r^^],^ and reorgani- 
zation energy. The maximum optimality and robustness for FMO is 
observed around the estimated value r^"^^ = Ins implying the sig- 
nificance of the time-scale separation between dissipation and trap- 
ping rates. We note that the ENAQT effect is ubiquitous at all rates 
of electron-hole recombination process. 



fied experimentally ilSTl . Up to this point, we have considered 
BChl 3 to be in the close proximity of RC by a trapping time- 
scale of about 1 ps. However, in this section we consider both 
of these parameters to be free, in order to explore the opti- 
mality and robustness of the ETE landscape as we vary the 
time-scale and geometrical constraints set by the RC trapping 
mechanism. 
In Fig. 
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we study the behavior of energy transfer effi- 
ciency landscape in various trapping time-scales and reorga- 
nization energies. It is evident that as the trapping rate be- 
comes very slow comparable to 100 ps or slower, the ETE 
drops significantly independent of the values of A. This can be 
understood intuitively as follows: the excitation has to wait on 
average so long for successful trapping to take place such that 
there will be an increasing chance of electron-hole recombina- 
tion as we are approaching time-scales comparable to exciton 
life time. If the trapping mechanism occurs within a time- 
scale of 1 ps, the ETE reaches to its expected maximum value 
of about 99%. Generally, one might expect that with increas- 
ingly faster trapping mechanisms the likelihood of dissipation 
to environment vanishes and energy transport approaches to 
the ideal case of having perfect efficiency. However, when 
the trapping rate becomes very fast on the order of lO^^ps or 
faster, the ETE also drops significantly, a result that might ap- 
pear counter-intuitive. In fact, overly rapid trapping leads to 
low efficiency via the quantum Zeno effect, as the rapid trap- 
ping effectively freezes the exciton dynamics and prevents it 
from entering the reaction center As a result, the finite exciton 
life-time eventually leads to complete dissipation of excitation 
to the environment in extreme limit of fast trapping of about 

The optimality and robustness of ETE versus both deco- 
herence and trapping rates using gradient and Hessian norms 
are presented in Fig. [TT](left panel). It can be observed that at 



A = 35cm~^ for the FMO, the ETE is optimal with a trapping 
rate of about 0.5ps. If the environmental interactions were 
stronger, a comparably faster trapping mechanism would be 
required to preserve such high level of efficiency. However, 
for small and intermediate system-bath interaction strength, 
where environment-assisted transport occurs, slower trapping 
rates become optimal, that is A and rtrap are not compet- 
ing processes anymore. This implies that a multi-parameter 
convergence of time-scales of the relevant physical processes 
might be required for light-harvesting complexes to operate 
optimally. We will address this issue in a broader context in 
the discussion section. From Fig. 10 it can be easily seen 
that ETE is very robust to variation of trapping rate at about 
1 ps time scale. In Fig. 11 we also illustrate the robustness 
with respect to both trapping and reorganization energy (right 
panel). For rather large values of A, there are certain regions 
that are highly robust to both parameters, but they are in fact 
suboptimal, as can be seen from noting their values in the left 
panel. 

To explore the optimality of ETE with respect to the lo- 
cation of reaction center, we consider the efficiency of other 
scenarios that the reaction center can be in the proximity of 



any other BChl sites. Figure 12 shows a cross-section of Fig. 
[2] in a fixed 7 — 50cm~^, for all possible trapping sites. To 
be unbiased with the respect to the initial state, we assume 
a maximally mixed initial state. It can be seen that the opti- 
mal curves belong to BChls 3 and 4 as expected, since they 
contribute highly to the lowest energy delocalized excitonic 
states. It should be noted that optimal environment-assisted 
quantum transport, and the two extreme regimes of quantum 
localization can be observed for all of these plots independent 
of the actual location of trapping. In other words, the behav- 
ior of the energy transport efficiency landscape and its depen- 
dence on a single governing parameter are not properties of a 
particular choice of trapping site in the FMO structure. 
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FIG. 14: Robustness of FMO transport efficiency witti respect to 
small variations of BChls locations, site energies and dipole orien- 
tations for 10000 samples: the Hamiltonian parameters are perturbed 
with site energies disorders of ±10cm~^, dipole-moment uncertain- 
ties of ±5°, and BChls spatial displacement of ±2.5A. The statis- 
tical distribution of 10^ random configurations shows a significant 
degree of robustness such that 99% of samples still preserve an effi- 
ciency of above 0.9. 



FIG. 15: Robustness of the FMO complex transport efficiency with 
respect to large variations in BChls site energies and dipole moment 
orientations for 10000 sample configurations: While the location of 
BChls are still slightly perturbed, similar to Fig. 
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of about ±2.5A, 

the dipole moments can take any arbitrary direction and site energy 
takes any value between zero and 500cm^^. This Histogram reveals 
that the relative distance of BChls is playing a crucial role in per- 
formance of these random light-harvesting complexes since 79% of 
them still hold ETE larger than 90%. 



For completeness, we also investigate the ETE landscape 
as a function of dissipation (loss) rate and reorganization en- 



ergy in Fig. 13 In our simulation of the FMO dynamics we 



have used the estimated value of r 



loss 



Ins. In Fig. 



however, we treat loss rate as a free parameter and we 
serve that for any stronger dissipation process, ETE would 
have a suboptimal and less stable behavior. Thus, even if all 
the other important parameters are within the optimal regime, 
a large time-scale separation between dissipation and trapping 
rate is still required to guarantee the highest performance for 
light-harvesting complexes. Fig. 13 also demonstrates that 



the existence of ENAQT is independent of a particular choice 
of dissipation rate. 



VI. OPTIMALITY AND ROBUSTNESS WITH RESPECT 
TO PARAMETERS OF FMO HAMILTONIAN 

So far, we have demonstrated that for the estimated Hamil- 
tonian elements of the FMO complex, the environmental pa- 
rameters and trapping rates are within the right set of val- 
ues leading to an optimal noise-assisted energy transfer ef- 
ficiency. Moreover, the performance of FMO is robust with 
respect to variations in such decohering and lossy processes 
and to uncertainties in initial conditions. However, in the con- 
text of exploring the role of coherent system evolution, it is 
not fully clear if the FMO internal Hamiltonian parameters 
have evolved to function optimally and fault tolerantly, de- 
spite disorders and thermal fluctuations. This issue has been 
examined for LHCII in Ref. 1 82 1 using semi-classical Pauli 
master equations to simulate the exciton dynamics. Here, we 
would like to use TC2 to explore how rare is the FMO geom- 
etry in terms of its efficiency, whether the specific spatial and 



dipole moment aiTangements of BChls are essential for such 
highly efficient functioning of this pigment-protein complex, 
and how robust these parameters are with respect to small and 
large perturbations in chromophoric distances, dipole moment 
orientations, and site energy fluctuations. Specifically, we ex- 
plore if the FMO closely packed structure plays any functional 
role, and illustrate a potentially important convergence of the 
relevant dynamical time-scales in the FMO energy transport. 
In the following section, we investigate the underlying physi- 
cal principle of quantum transport in more generic multichro- 
mophoric structures beyond the FMO geometry. We examine 
the efficiency of small light-harvesting complexes as functions 
of the compactness level, number of sites, reorganization en- 
ergy, exciton energy gaps, and multi chromophoric connectiv- 
ity. 

The Frenkel exciton Hamiltonian for a multichromophoric 
system is expressed as: 



Hs = }_^e,\j}{j\ + J,,k\j}{k\, 



(10) 



in which Jjk are Coulomb couplings of the transition densities 
of the chromophores. 



jk 



(/^j ■ R-ifc)(Mfc • Rjfc)), (11) 



where TLjk denotes the distance between site j and k, and fij 
is the transition dipole moment of chromophore j fSl. We first 
study the robustness of free Hamiltonian parameters within 
the proximity of the estimated values for FMO as given in the 
appendix A. 
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A. Robustness of FMO Hamiltonian 



Figures 14 and 15 demonstrate the robustness of the FMO 
structure to variations in its internal parameters. Figure [l4|il- 
lustrates that FMO efficiency does not drop drastically with 
respect to variations in the dipoles orientations, site energies, 
and Bchls distances close to the neighborhood of the estimated 
values. Specifically, from 10000 random samples of FMO 
with spatial uncertainty around each Bchl location of about 
±2.5 A, dipole moments orientations variations of ±5°, and 
site energy static disorder ±10 cm^^, 97% of configurations 
have efficiency of 95% or higher. This demonstrates a sig- 
nificant degree of robustness with small perturbations. In or- 
der to separate the influence of spatial coordinates from angu- 
lar dipole orientations and disorders, we allow the latter two 
parameters to take arbitrary values from a large range while 
keeping Bchl locations uncertainties to be limited by ±2.5 



A. We observe in Fig. 15 that ETE remains relatively robust 



with 79% of random 10000 configurations have still efficiency 
of 90% or higher. This is rather counterintuitive considering 
huge freedom that we have accommodated in the dipole mo- 
ment arrangements and site energies. These results indicate 
that spatial degrees of freedom is a dominating geometrical 
ingredient of the FMO structure and might play a key physi- 
cal role in its performance. 




Reorganization Energy (cm ) 



Scaling Factor k 



FIG. 16: Dependency of ETE on compactness level of the FMO com- 
plex. The FMO chromophoric spatial structure is scaled with a factor 
between 0.5 and 5. The ETE manifold is plotted as a function of the 
scaling factor and reorganization energy. For all levels of compact- 
ness, A = 30 — 40cm~^ yields the highest ETE. A more compact 
complex shows a higher degree of robustness with the respect to vari- 
ation in reorganization energy. The scaling factor 1 corresponds to 
the FMO complex leading to a nearly optimal energy transfer effi- 
ciency. 



B. Searching for FMO-type Hamiltonians 

To explore the role of the free Hamiltonian parameters, we 
first investigated energy transfer efficiencies of 200,000 ran- 
dom 7-site chromophoric configurations embedded within a 
sphere of diameter 200A, with random initial and target sites, 
and with environmental parameters the same as the estimated 
values of FMO |64|. We bounded the nearest distance among 
chromophores to be larger than 5A, since the dipole-dipole 
approximation breaks down below this limit [il8|. The dipole 
moment of each chromophore can take any orientation and the 
site energy of each chromophore can take any value between 
zero and 500cto~^ chosen from a random distribution. We ob- 
served that only 10% of these configurations have efficiency 
larger than 50% and only 0.1% have efficiency comparable 
to FMO. Thus, one might conclude that the FMO structure 
is one of the rare geometries that can support this extremely 
high efficiency. However, as we argue below this conclusion 
is inaccurate. 

We searched over various geometrical pattern correlations 
among the top 0.1% high-end efficient random structures to 
reveal the underlying principle(s) that could give rise to their 
high efficiency. We found that the distance between initial and 
target sites in most of these structures vary between 6 — lOA. 
Thus most of these complexes are rather trivial optimal solu- 
tions that are not particulary interesting for quantum transport 
studies. Indeed each of those conformations essentially acts 
as a dimer including a donor and an acceptor chromophore 
that are strongly coupled via dipole-dipole interaction in the 
near field. This approach for optimality is obviously not us- 
able for quantum transport in large natural or engineered light- 



harvesting antenna or macromolecular systems with typically 
significant spatial separation of initial excitation and trap- 
ping sites. In order to distinguish potential non-trivial opti- 
mal solutions hidden in these results; i.e., those involve quan- 
tum/classical jumps in a multi-path molecular networks, we 
consider the distance of donor and acceptor to be bounded by 
the size of the complex for the rest of our studies in this paper 
Our results presented in Figs. 14 and 15 clearly indicate 



that the relative Bchl locations play a major role in the overall 
performance of the FMO complex. Thus a potentially sig- 
nificant parameter of relevance is the compactness of a given 
pigment-protein complex. We further explore this feature by 
introducing a single compactness parameter by rescaling the 
Bchl distances by a factor k. We plot ETE as a function of 
compactness level k varying by an order of magnitude from 
0.5 to 5. To explore any potential interplay of environmen- 
tal interactions with this particular internal degree of freedom, 
we also simulate this size dependent supersssion/enhancement 
of ETE in various reorganization energy. A, in Fig. 16 It 
can be seen that transport efficiency drops significantly by ex- 
panding the FMO structure and has slightly enhanced value 
for contraction ratio of 0.5. As above, in the discussion of 
the evolutionary tuning of reorganization energies and envi- 
ronmental correlation times, it appears that nature has taken a 
rather minimal approach in designing FMO. By utilizing the 
right density level of chromophoric structure and strength of 
environmental interactions, the evolutionary process has cre- 
ated a nearly optimal energy transfer channel, but that process 
apparently avoided the extra work necessary for further com- 
pression to achieve just 1% or less enhancement in ETE. By 



inspection of the ETE variation with resect to A in Fig. 16 



one observes that the ENAQT phenomenon is scale invariant 
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for the FMO-like structure. 

The above analysis reveals an important effective physi- 
cal parameter; however, one cannot rule out the possibility 
of other underlying dynamical processes and/or geometrical 
correlations that could also lead to higher ETE in very dif- 
ferent compactness regimes. In other words, the compactness 
appears to be a sufficient condition for high ETE in multi- 
chromophoric complexes but most hkely it is not a necessary 
condition. 



C. Convergence of time-scales for FMO 

We illustrate a summary of our main results on the FMO 
complex energy transport dynamics in the non-Markovian and 



35(K) ^ . 66.7(fs) 



Bath temporal 
correlations i , 



non-perturbative regimes in Fig. 17 This figure indicates that 
an effective convergence of time-scales for various dynamical 
processes exists at the intermediate and optimal regimes (de- 
noted by a white circle) in the energy transfer landscape rep- 
resented by a circular surface. The radius of this circular area 
illustrates the reorganization energy in logarithmic scale in the 
unit of time from 33.3ps (1 cmT^) to 66.7/s (500cto~^). The 
circular area consists of four different regions in four quarters, 
the angular degrees of freedom in each quarter is associated 
to a different physical quantity. The first quarter represents 
the degree of the memory in the bath or the non-Markoviaity 
of the environment in logarithmic scale, that is equivalent to 
bath correlation time, varying from about Q.&Ips to 66.7/s. 
The second quarter illustrates variations in the bath tempera- 
ture from 35 to 350° iiT. The third region shows six order of 
magnitude trapping time-scales from Ins to Ips in a logarith- 
mic scale. The last quarter demonstrates compactness level 
of the Frenkel exciton Hamiltonian of the FMO rescaled by 
a factor from 0.5 to 5. The estimated parameter of the FMO 
is denoted by a white bar in each region. In all quarters the 
efficient neighborhood of ETE are somewhere in the medium 
range, denoted by a white circle at about 1 ps time-scale cor- 
responds to reorganization energy of the FMO. The optimal 
areas are far from extreme low and high reorganization ener- 
gies in the center and perimeter of the black circle respectively 
where quantum localization effects can substantially suppress 
ETE. The key point of this figure is that all the relevant use- 
ful transport time-scales converge just around Ips, far from 
the dissipation time-scale of Ins. We will discuss the impli- 
cations of these results in a broader context of optimality and 
robustness of complex quantum systems in the next section. 



VII. DISCUSSIONS ON THE CONVERGENCE OF 
TIMESCALES: QUANTUM GOLDILOCKS PRINCIPLE 

Our numerically study indicates optimality and robustness 
of FMO with respect to a variety of parameters. For FMO, 
in all cases, we find that the estimated operating values sit 
nearly at the optimal point for energy transport efficiency, 
and that there is a relatively wide range of parameters about 
that point for which FMO is close to optimal. Moreover, we 
note that the estimated time/energy scales associated with this 




FIG. 17: A summary of our main results on the optimality and robust- 
ness of FMO energy transfer dynamics illustrating a convergence of 
time scale around Ips. The radial axes indicate the reorganization en- 
ergy in logarithmic scale which is essentially the main environmental 
parameter quantifying the strength of system-bath interactions. An- 
gular coordinate in each quarter represents a different physical quan- 
tity. The first three quarters are associated to three other important 
environmental parameters including bath correlation time-scale, tem- 
perature, and trapping time-scale. The last region highlights the role 
of compactness as the most significant internal parameter of FMO 
structure quantified by a rescaling factor in the chromophoric rela- 
tive distances. To improve the clarity the color bar is rescaled in each 
quarter. The estimated values of the FMO complex coincide in the 
intermediate values of both angular and radial coordinates, denoted 
by white lines (around Ips), leading to optimal and robust ETE. 



point - couplings between sites, reorganization energies, de- 
coherence times, environmental correlation time, and trapping 
time-scale- all lie in the same intermediate range, on the or- 
der of a picosecond. Nevertheless it is instrumental for the 
FMO performance that life-time of exciton has a sufficiently 
longer time-scale. For alternative configurations, beyond the 
specific FMO model, we observe that their chromophoric den- 
sity, geometrical arrangements, and coupling strengths need to 
be in tune with environmental interactions to lead to optimal 
performance. Here, we briefly discuss a design principle that 
fundamentally relies on this convergence of time scales. 

A common design technique in engineering is to keep the 
desired system as simple as possible. The rationale is that 
overly complex designs are typically not robust: the more 
complex the design, the harder it is to predict its behavior, 
and the more pieces there are to fail. Now, FMO and other 
photosynthetic complexes are definitely not simple: they con- 
sist of multiple parts, and, as our studies show, enlist a wide 
variety of quantum effects, tuning those effects to attain effi- 
cient and robust operation. Apparently, these complexes do 
not obey the fundamental engineering philosophy of keep it 
simple. There is, however, a more sophisticated version of 
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the design principle in system engineering, as overly simple 
designs may be robust, but they typically also fail to attain 
their desired functionality. This advanced version of the de- 
sign principle declares that to the extent that it is possible to 
add more features to a design to enhance functionality with- 
out reducing robustness, it is desirable to do so. That is, given 
available technologies, there is typically a design that attains 
a level of complexity that is 'just right' - not so simple that it 
fails to attain its goal, but not so complex that it loses robust- 
nesses. The design strategy that aims for the just-right level of 
complexity is sometimes called the Goldilocks principle for 
complex systems |[83l . Let's look at how this principle ap- 
plies in the case of the structure of FMO, a complex quantum 
system that has been evolved by natural selection. 

The Goldilocks principle for complex quantum systems - 
the quantum Goldilocks principle- states that to improve effi- 
ciency, quantum effects could be used, as long as by doing so 
it does not make the operation of the system more fragile. Our 
simulations show that this is the case with FMO: a wide vari- 
ety of parameters have been tuned to their optimal value, and 
that optimal value is robust against variations in those param- 
eters. The quantum Goldilocks principle essentially explains 
why it is beneficial for certain parameters that involve time or 
energy to converge to the same time/energy scale. 

When the time scale for one effect is very different from 
the time scale for another, then the first effect can typically be 
regarded as a perturbation on the other. For example, if the 
environmental interactions are much slower than the system 
internal dynamics, and/or if the correlation time of the envi- 
ronment is much shorter than the coherent tunneling rate of 
an exciton from chromophore to chromophore, then the effect 
of the environment can be treated in the perturbative and/or 
Markovian limits. If one effect is to interact strongly with an- 
other effect, by contrast, then the time/energy scales of the 
two effects should be similar. As we have shown here, the 
maximum efficiency of quantum transport is attained when the 
key decoherence parameter AT/7 is on the same order as the 
free Hamiltonian coherent time-scale g. As shown elsewhere 
Il29l[3ni33ll67l environmentally assisted quantum transport is 
an important effect. This convergence of time scales is a nui- 
sance for quantum simulation, as perturbative methods break 
down. While this phenomenon makes the efficient simula- 
tion of such systems very difficult, the convergence of time 
scales is quite useful for nature, however, as the two effects 
can now interact strongly with each other to produce a signif- 
icant enhancement in efficiency. Of course, too strong inter- 
action could also produce a significant decrease in efficiency. 
But for systems undergoing natural selection, one typically 
expect that the strong interaction is tuned at the right level to 
give rise to a beneficial effect. 

The results of this manuscript illustrates the quantum 
Goldilocks principle manifestation in the FMO complex 
and other small-size mutlichromophoric configurations: the 
time/energy scales have converged within the range where the 
different quantum effects and structural parameters at work 



can interact strongly with each other to give rise to possi- 
bility of high efficient energy transport. Moreover, this va- 
riety of quantum effects can be combined to attain high ef- 
ficiency without sacrificing robustness. The convergence of 
time scales makes these complexes difficult to model, but 
highly efficient in certain instances. We anticipate that a simi- 
lar convergence of time/energy scales could appear in larger 
biomolecular complexes and in biological processes where 
quantum mechanics might play an important role. 



VIII. CONCLUSION 

In this work, we studied the structural and dynamical design 
principles of excitonic energy transfer in the Fenna-Mattews- 
Olson (FMO) complex. Our numerical simulations demon- 
strate that the natural structure of FMO pigment-protein com- 
plex lead to a highly efficient and robust photosynthetic en- 
ergy transfer wire. We characterized the ETE landscape by 
three main regions: weal localization, environment-assisted 
quantum transport, and strong localization, and identified the 
scalar A = ^ as the key parameter to cross between these 
regions as one hikes over the landscape. 

We showed optimality and robustness of energy transport 
efficiency in FMO with respect to variations in all the main ex- 
ternal parameters dictating the dynamics. In particular, we ex- 
plored the protein-solvent environment factors: system-bath 
coupling strength, bath memory, temperature, and spatial cor- 
relations and dissipation; the light-harvesting antenna factor: 
electronic state initialization; and reaction center factors: ex- 
citon trapping rate and location. Next, we investigated the 
performance and sensitivity of the FMO with small and large 
variations in its internal parameters: chromophoric spatial lo- 
cations, dipole moment orientations, site energies, and chro- 
mophoric density. The exciton trapping process in FMO com- 
plex shows to be significantly robust with respect to perturba- 
tions in dipole moment orientations and site energies due to 
its compact structure. 

We believe that a key design principle for achieving op- 
timal and robust quantum transport performance is to allow 
for the convergence of energy/time scales for the contribut- 
ing internal and external parameters. This phenomenon can 
be regarded as an examples of a Goldilocks principle in the 
quantum regime. 
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FIG. 18: The top view of FMO energy transfer landscape in the pres- 
ence of an Ohmic model of bath spectral density as a function of re- 
organization energy and bath frequency cutoff at room temperature. 
The ETE predicted by this model is below the estimate values by a 
Lorentzian model and therefore it cannot explain near ideal perfor- 
mance of FMO complex in exciton transport. Nevertheless, here we 
can also observe the three distinct regimes of weak and strong local- 
izations close to 7 and A axes, and the intermediate optimal ENAQT, 
separated by straight lines. This plot confirms the role of parameter 
AT/7 &s ™ effective governing parameter. 
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In this work we use the following free Hamiltonian for the 
FMO complex, given in Ref. [4|: 
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The estimated values of dipole moment orientations and po- 
sitions of the Mg atoms, representing the location of Bchls, are 
extracted from the pdb file of the FMO complex |84|. These 
data can be summarized in the following table: 

The Bchl-Bchl coupling in FMO is dipole-dipole interac- 



FIG. 19: Top view of the ETE landscape at 11° K for a bath with a 
Lorentzian spectral density. The general pattern of an optimal trans- 
port region in between two limits of localized low-efficient transport 
regimes, divided by rather straight lines, is observed here as well. 
However, due to low-temperature effects captured in the parameter 
AT/7, the adversarial high reorganization and slow bath effects be- 
come less pronounced leading to a higher ETE than those in Fig. [3] 
in the corresponding limits. 



tion 

C 3 
J]k = ^(Mj ■ Mfe - ^(/^j ■ R.jfe)(Mfe • Rjfc)), (Al) 

for which we choose the constant C|/ip = 134000 cm~^ 

ma. 
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TABLE I: Spatial location of Bchls and their dipole moment orienta- 
tion. 



Bchl 


x(l) 


y(A) 


z(A) 


e 


<!> 


1 


28.032 


163.534 


94.400 


0.3816 


-0.6423-H7r 


2 


17.140 


168.057 


100.162 


0.067 


0.5209-^7r 


3 


5.409 


180.553 


97.621 


0.1399 


1.3616-^7r 


4 


9.062 


187.635 


89.474 


0.257 


-0.6098-H7r 


5 


21.823 


185.260 


84.721 


-0.1606 


0.6899-^7r 


6 


23.815 


173.888 


82.810 


-0.4214 


-1.4686-H7r 


7 


12.735 


174.887 


89.044 


0.578 


-1.0076-^7r 



Appendix B: Energy transfer efficiency for Ohmic bath and 
cryogenic temperature 

Here, we demonstrate that separation of ETE landscape, as 
a function of the parameter AT/7, into the various regions 
with distinct quantum transport efficiencies is not a mere prop- 



erty of either Lorentzian bath or high temperature limit. A 
top view of FMO energy transfer landscape is shown for an 
Ohmic spectral density at 298° X and Lorentzian spectral den- 



sity at 17° K in Fig. 18 and Fig. 19 respectively. In each fig- 
ure, one can distinguish three different regions including low- 
efficient weakly localized limit, optimal EANQT, and low- 
efficient strongly localized limit, similar to the results pre- 
sented in Fig.|3] 

It should be noted that a bath with an Ohmic regularized 
spectral density, J(a;) = X{lli/^) exp(— aj/7), has often been 
employed in modeling the effect of bath fluctuations on the 
spectroscopic readout of an FMO sample ||4]|72l. Figure 18 
presents the top view ETE landscape as a function of bath re- 
organization energy and cutoff frequency with an Ohmic den- 
sity. The observed high efficient region stays lower than its 
Lorentzian counterpart in Fig. |3] therefore is unable to explain 
the high efficiency of the FMO complex. This can be seen 
as a confirmation for the theoretical modeling of the solvent- 
protein envkonment with Lorentzian spectral densities 1.57 J . 



